MULTI-POINT TAYLOR EXPANSIONS OF ANALYTIC 

FUNCTIONS 



JOSE L. LOPEZ AND NICO M. TEMME 



Abstract. Taylor expansions of analytic functions are considered with re- 
spect to several points, allowing confluence of any of them. Cauchy-type for- 
mulas are given for coefficients and remainders in the expansions, and the 
regions of convergence are indicated. It is explained how these expansions can 
be used in deriving uniform asymptotic expansions of integrals. The method is 
also used for obtaining Laurent expansions in several points as well as Taylor- 
Laurent expansions. 



In deriving uniform asymptotic expansions of a certain class of integrals one 
encounters the problem of expanding a function, that is analytic in some domain 
Q of the complex plane, in several points. The first mention of the use of such 
expansions in asymptotics is given in pQ, where Airy- type expansions are derived 
for integrals having two nearby (or coalescing) saddle points. This reference does 
not give further details about two-point Taylor expansions, because the coefficients 
in the Airy-type asymptotic expansion are derived in a different way. Other men- 
tions of the use of such expansions in asymptotics is given in and [H]- In jjj, 
two-point Taylor expansions are used with applications to Airy-type expansions of 
parabolic cylinder functions. In we used two-point Taylor expansions to derive 
convergent expansions of Charlier, Laguerre and Jacobi polynomials in terms of 
Gamma, Hermite and Chebyshev polynomials respectively. 

To demonstrate an application in asymptotics of multi- Taylor expansions we 
may consider contour integrals of the form 



where a is a vector of parameters, a — (ai, . . . , a-i) and the phase function f(z, a) 
has m saddle points Z\, Z2,...,z m . The asymptotic behaviour of these integrals for 
large values of A is determined by the saddle-point structure of the phase function 
chap. 7, sec. 6]. One method for obtaining an asymptotic expansion of this 
integral for large values of A is based on expanding g(z) at the saddle points of the 
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phase function, 

oo 

(1.2) g(z) = + atz + . . . + a m ^z m -%z - Zl ) n (z - z 2 ) n ■ ■ ■ (z - z m ) n 

and substitute this expansion into (|l.lfl . When interchanging summation and in- 
tegration, the result is a formal expansion in m series in terms of functions related 
with the functions 
(1.3) 

F n . k (X;a) = f z k (z- Zl ) n (z-z 2 ) n ■■■(z-z m ) n e- x ^' a Uz, k = 1, 2, . . . , m-1. 

In 0, these functions F n k(\;a) are the Airy functions, whereas in [5] these func- 
tions are the Gamma function, or the Hermite or Chebyshev polynomials. 

In a future paper we will use multi-point Taylor expansions in the asymptotic 
analysis of integrals arising in diffraction theory, such as the Bessel function integral 
(see and U) 

(1.4) J(x,y)= / U (yt)e^ t4+xt2 Ut, 

Jo 

which is related to the Pearcey function 

-oo 

The Taylor-Laurent expansions will be used to study integrals with two saddle 
points and a pole of the integrand. Other applications in asymptotics include the 
study of Hermite-Pade approximations to the exponential function; in [2] integrals 
are considered with three saddle points. 

In a recent paper |S] we have introduced the theory of two-point Taylor ex- 
pansions, two-point Laurent expansions and two-point Taylor-Laurent expansions. 
The purpose of the present paper is to generalize that theory from 2 to m points, 
m > 2. We give details on the region of convergence and on representations in terms 
of Cauchy-type integrals of the coefficients and the remainders of the expansions. 
Earlier information on this type of expansions is given in [S], Chapters 3 and 8. 
The theory of several-point Taylor expansions was already formulated in Chapter 
3 of Walsh's book, although in a different setting. Chapter 8 of presents also a 
theory of rational approximation of analytic functions, but is different from the the- 
ory of multi-point Laurent and Taylor-Laurent expansions presented here. Whereas 
the multi-point polynomial approximation of Chapter 3 may be reformulated as a 
multi-point Taylor approximation, the rational approximation of Chapter 8 can not 
be written as a multi-point Laurent or Taylor-Laurent approximation. For more 
details, see Section 5. 

2. Multy-point Taylor expansions 

We consider the Taylor expansion of an analytic function f(z) in several points and 
give information on the coefficients and the remainder in the expansion. In what 
follows empty sums and derivatives of negative order must be understood as zero 
and empty products as one. We will deal with the following set of points: 

Definition 2.1. We define the set 

(2.1) S = {zi,Zl, . . . , Z\\ Z2, 22, ■ ■ • , 22, ■ ■ • 5 Zp, Z p , . . . , Z p } 
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of to points consisting on p different points zi, Z2,..., z p (zi ^ Zj if i j), each Zj 
repeated rrij times: mi + to,2 + . . . + m v = to. 

For clearity in the exposition, we first introduce the multi-point Taylor expansion 
for to different points Z\, Z2, ... z m (m = p, m.j = f ) in Theorem 12.31 In Theorem 
12. 51 we assume that the points z\, z%, . . . z m may coalesce. We will need the following 
elementary lemma. 

Lemma 2.2. Given z, w, £ (D, take m different points z\, Z2,..., z m in (D and 
define 

(2.2) H m ( W , z m ) . nLi^-^-nLi^-^) . 

w — z 

Then 

(2.3) H m (w, z;zi,...,z m )='J2 ^ fc ~ 1 ' M3 i r m -7T~~T- — • 

Proof. The numerator of H m (w, z; z±, . . . , z m ) is a polynomial of degree to in the 
variable w that vanishes at w = z. Therefore, H m (w, z; z\, . . . , z m ) is a polynomial 
of degree m — I in the variable w. Let P m (w 7 z; zi, . . . , z m ) denote the function 
at the right-hand side of l|2.3(l . which is also a polynomial of degree to — I in the 
variable w. Moreover, 

m 

(2.4) H m (z s ,z;zi, ... ,z m ) = J][ (z - Zk) = P m (z s ,z;zi, . . . , z m ) 

k=l,k=jLs 

for s — 1, 2,. . . to. Hence, 

(2.5) H m (w, (u>, z; Zi, . . . , z m ). 

□ 

Theorem 2.3. Let f(z) be an analytic function on an open set f2 C C and S C f2 
with S consisting on m different points (to = p). Then. f(z) has the multi-point 
Taylor expansion 

N-l m 

(2.6) f(z) = J2 ~ z 0" + rjv(«), 

n=0 fe=l 

where q nt m{z) is the polynomial of degree m — 1 



(2.7) q n ,m{z) = y^«n, 



nfc=i,*#j(«-«fc) 



j =1 llk=l,kjtj\ Z j _ 

and i/ie coefficients a n j of this polynomial are given by the Cauchy integral 

OR\ n =— f /H dw 

( j ^~2mJ c ( W - Zj m: =1 ^-^r 

The remainder term r^(z) is given by the Cauchy integral 



(2.9) r N (z) __ 

Zm 



.N 

1 f f{ w ) dw 



c (w-z) [lTLi(w ~ z k )f 



]J(Z-Zk) 
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The contour of integration C is a simple closed loop which encircles the points z\, 
Z2,.-.,z m (fora n j) andz, z\, Z2,...,z m (for rj\r(z)) in the counterclockwise direction 
and is contained in fl (see Figure 1 (a)). 

The expansion \2. b)) is convergent for z G O rn , where: 
(2.10) 

rn ( m \ 

O m { : <>. ! [ : :/, • r\. r = ltif,,. ( < J J w :/. > • 

k=l lfc=l J 

That is, y2,.b)) is convergent for z inside the lemniscate JlfeLi \ z ~ z k\ = r (see Figure 
2; if m = 1 this domain is a disk; if m = 2 this domain is bounded by a Cassini 
oval). In particular, if f(z) is an entire function (fl = (B), then the expansion \2.(>\) 
converges V z G (27. 



(2.H) f(z) 



Proof. By Cauchy's theorem, 

1 f f(w) dw 
2ni J c w — z ' 

where C is the contour defined above (Figure 1 (a)). We write 

) 1 

w - z \ [ k=1 (w - Zk) 1 - u 
where H m {w, z; Z\, . . . , z m ) is given in 12.211 and 

(2.13) u ^f- Z "\ . 

Now we use Lemma l2~2l and introduce the right hand side of (|2.3|l and the expansion 

(2.14) -J- = J2 « B + J— 

1 — u i 1 — u 

n=0 

in (|2.12|) and this in (|2.11|) . After straightforward calculations we obtain formulas 
(l2~Tni-(l2~31). 

For any z € O m , we can take a contour C in f2 such that 

m m 

(2.15) J] \z-z k \ < Y[ \w-z k \, VioeC 
fe=i fe=i 

(see Figure 1 (b)). On this contour is bounded by some constant C: \f(w)\ < 

C. Introducing these two bounds in (|2.9|l we see that limjv->oo tn(z) = and the 
proof follows. □ 

We need the following lemma to consider the case of coalescing points in the 
set S. 

Lemma 2.4. Given z,w G (B, take m different points z\, z%, z m in C, all 
different from w too. Let those m points to coalesce at z m , say. Then 

(2-16) hm > —jjpl = > 

z 1 ,z 2 ,..,z m _ 1 ^z m ^— ' ( W - Zi) L_1 u.Ui - Z k ) *7~L 



3- 
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FIGURE 1. The case m — 3. (a) Contour C in the integrals 12. 8f) 
and (|2.9I) . (b) For z £ O nll we can take a contour C in f2 which 
contains O m inside and therefore, JlfeLi \ z ~ z k\ < IlfcLi \ w ~ z k\ V 
w £ (D. 







r = n 



Figure 2. Shape of the "lemniscate domain" O m for to = 3. It 
depends on the size of the parameter r defined in H2.10(l . In these 
pictures |z 2 — z 3 | < \z\ — z 3 |, \z\ — z 2 | and r a > r^> r c . 



Proof. We first note that the identity 

n j— 1 n 

(2.17) En^-^) n («i-*i)=o 

j=i i=i i=i+i 
holds for any set of points zi, z 2 ,..., z„, n > 1. It may be cheked in the following 
way: we take the two first terms of the sum, which gives 

(2.18) {z n - z 2 )(zi - z 3 )(zi - za) ■ ■ ■ (zi - z n -i)(zi - z n ). 
Next we add to this the third term of the sum, which gives 

(2.19) (z„ - Z 2 )(z n - Z 3 )(zi - Z 4 ) ■ • ■ (zi - Z n -l)(zi - Z n ). 

We continue this process untill we add the n — 1-th term of the sum, obtaining 

(2.20) (z„ - z 2 )(z„ - z 3 )(z n - Zi)- - ■ (z„ - z„_i)(zi - z„). 

But this is just the last term of the sum with opposite sign. 
Using the above identity we have 



fc-i 



(2.21) Y{{ z s-zi) 



i=i 



3-1 



1 [ (z k -zi) = 



l=s+l 



j=k l=k 1=3 + 1 

for any s = 1, 2, 3, to and any k = 1, 2, 3, m with k ^ s. Then 

s j — 1 m m J — 1 m 

(2.22) e - *o n (** - = e - n (** - *o = ° 

for any s,k = 1, 2, 3, m with k =/= s. 
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Now, for every s = 1,2,3, ...,m, we define the following polynomials of degree 
to — 1 in the variable z: 

m m j—1 m 

(2.23) R s {z)= J] (z- Zl ), S s {z) =Y,n(z ~ zt) [] (z a -zt). 

l=l,l=£s j = l 1 = 1 l=j + l 

The zeros of R s (z) are Zk for k — 1,2,3, ...,to, k ^ s and from (|2.22() . S s (zk) = 
for k = 1,2, 3, to, k ^ s. Moreover, the leading coeficient of R s (z) and of S s (z) 
coincide. Therefore, R s (z) = S s (z) for s — 1,2, 3,..., to. 

Finally, define the following polynomials of degree to — 1 in the variable w: 

(2.24, p(,..)-£ n^u-^)rE..^,(— ».) | 

m J—1 m 

(2.25) Q( W ,z) = ^]J(z-z fc ) [J (ui-* fc ). 

For every s = 1,2,3, to we have P(z s ,z) = R s (z) and Q(z Sl z) ~ S s (z). But 
R s (z) = S s (z) and therefore P(w, z) — Q(w, z). Then, 

"h( w - z j) nr=i,fc#j ( z i ~ zfc ) nr=i («> - **o 

(2-26) m 

Taking the limits z\, z%,.., z m _i — » z m in the left- and right-hand sides of these 
equalities, we obtain the desired result. □ 

Theorem 2.5. Let f(z) be an analytic function on an open set Q C (D and Scfl. 

Then, f(z) has the multi-point Taylor expansion 

N-i p 

(2.27) f(z) = J2 qn,m(z)T[(z ~ Zk) nmh +r N (z), 

n=0 fe=l 

where q n ,m{z) is the polynomial of degree m — 1 

(2.28) q ntm (z) = Y[V~ M (? - z,\^ ^ a " J '' (z ~ ^ 

j=1 llk=l,k^j\ Z J Z k) l=Q 

and the coefficients a n jj of this polynomial are given by the Cauchy integral 
, 9 on\ - 1 f /M dw 

The remainder term rjv(z) is given by the Cauchy integral 



(2.30) r N (z) = 7T~- f 
2tti J c 



f(w) dw 



c{w~z)[Y\l =1 {w-z k r*] N 



n (*-**)' 



N 



Lfe=i 

The contour of integration C is a simple closed loop which encircles the points z\, 
z%,...,z v (for a n ,j,i) and z, z\, z%,...,z p (for Tn{z)) in the counterclockwise direction 
and is contained in Q (see Figure 1 (a)). 
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The expansion \2.2T[ is convergent for z 6 O v : 
(2.31) 

v 

O p = {zen, Y[\z-z k \ m * <r}, r 
k=i 



Inf 



n 



w - z k \ 



that is, inside the lemniscate n?_i \z — Zfe|" lfc = r. In particular, if f(z) is an entire 
function (£1 = (B), then the expansion converges V z £ (D. 

Proof. If all the points in S are different, we have from (|2.7(l and (|2.8(l 

nr=i,fc#j( z - z fc) 



(2.32) 



1 

2vr7 



f(w)dw__ x - 



This last sum may be also decomposed in the form 
(2.33) 



E 



E 



E 



]=m p - 



lYk=l.k^( Z ~ Z k) 

1+1 ( w - z j)rifcLi,fe^(^ -«*)' 



Now let the first mi points to coalesce to z\, the second mi points to coalesce to Z2, 
and so on, and apply Lemma 12^31 to every one of the p sums above to obtain i|2.27[l . 
(12.28(1 and (|2.29|l . Equation l|2.30() follows from (|2.9I) . The proof of the convergence 
of f2~2"?jl in the region O p is a straightforward generalization of the correponding 
proof in Theorem 12.31 □ 

2.1. Explicit forms of the coefficients. Formula (|2.29j) is not appropriate for 
numerical computations. A more practical formula to compute the coefficients of 
the above multi-point Taylor expansion is given in the following proposition. First 
we have a definition: 

Definition 2.6. Let f(w) be analytic at w; then for n — 0, 1, 2, . . . the differenten- 
tial operator D™f(w) is defined by 



(2.34) 



1 d n 
n! dw 1 



7H- 



Proposition 2.7. The coefficients a n ji, for n = 1,2,3,. 



J 



1,2,. ..,p, I 



0, 1, . . . , raj — 1 in the expansion 



are also given by the formula: 



jjnrrij 



(2.35) 



+ 



E D T k ~ x 



Proof. We deform the contour of integration C in equation H2.29fl to any contour 
of the form C\ U C% U • • • U C p , also contained in Q, where C k , k — 1,2, ... ,p, is 
a simple closed loop which encircles the point z k in the counterclockwise direction 
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(a) (b) (c) 

Figure 3. Integration contours Ck for p = 3 and q = 1. 
(a) The function Hs=i s^k( w ~ z s )~ nms f{w) is analytic in- 
side Cfc for k = 1,2,..., p. (b) The functions Ils=i s ^k( w ~ 
z s )~ nm * 9k{w) and EKUi.^^ ~ •z s ) ( ™ +1)ms .9fe(w) are analytic in- 
side Cfe for fc = 1,2,..., p. (c) The functions YV s =i s ^k^ w ~ 

z s )- nm °9k(w), nLw*(«'-*.)" nm ' rE= g +i,^*(«'-^) Am 'flfcW 

andnLi >? ^(^-^)- (n+1)ms x m= g+ i, s ^(^-^) ( " +1)ms 3fcH 
are analytic inside Cfc for k = 1, 2, . . . ,p. 



and does not contain any other point Zj, j = 1,2, ... ,p, j ^ k inside (see Figure 3 
(a)). Then, 

1 f /H dw 

1 /■ 



from which equation (|2.35|) follows. □ 



2.2. Multi-point Taylor polynomials. In Theorem 12.31 we have assumed that 
the function f(z) is analytic in fl. If f(z) is not analytic in f2 but has a finite number 
of derivatives at z\, zi,...,z p , we can still define the multi-point Taylor polynomial 
of the function f(z) at z\, Z2,...,z p in the following way: 

Definition 2.8. Let z be a real or complex variable. If f(z) is Nm,k — 1— times 
differentiable at z\, Z2,...,z p , we define the multi-point Taylor polynomial of f(z) at 
the points of S and degree mN — 1 as 

N-l p 

(2.37) P N (z) = J2 1n,m(z) j[(z z k ) nm \ 

n=0 fe=l 

where q n ,m{z) is the polynomial of degree m — 1 

p U p (7 - 7, ) mk mj ~ 1 

(2.38) qn , m (z) = £ W(z - E " 

J = 1 llk=l,k^ 3 \ Z 3 Z k) l=Q 

and the coefficients a nji i are given in Ij2.35|l . 

Proposition 2.9. In the conditions of the above definition, define the remainder 
of the approximation of f(z) by Pn{z) at the points of S as 

(2.39) r N {z) = f{z)-P N {z). 
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Then, (i) rjv(z) = o(z — Zk) Nmk 1 as z — > Zk, k = 1,2,..., p. (ii) If f(z) is 
Nmk — times differentiable at Zk for some k, then tat(z) = 0(z — Zk) Nnik as z — > Zk- 

Proof. The proof is trivial if f{z) is analytic at every z\, z 2 ,...,z p by using i|2.3U[l . 
In any case, for real or complex variable, the proof follows by using l'Hopital's rule 
and □ 

Remark 2.10. Observe that the Taylor polynomial of f{z) at the points of S and 
degree mN — 1 is the Hermite's interpolation polynomial of f(z) at Z\, z%,...,z p with 
data f'(z k ),...,f( Nm "- 1 \z k ),k = l,2,..., P . 

3. Multi-point Laurent expansions 

In the standard theory for Taylor and Laurent expansions much analogy exists 
between the two types of expansions. For multi-point expansions, we have a similar 
resemblance in the representations of coefficients and remainders. 

Theorem 3.1. Let Qq and Q be closed and open sets, respectively, of the complex 
plane, and Qq C C (B. Let f(z) be an analytic function on D, \ Oq and Z\, 
Z2,...,z p 6 fio (That is, S G Qo)- Then, for any z e f2\f2oj f(z) has the multi-point 
Laurent expansion 

N-l p N-l p 

(3.1) }{z) = J2 Qn,m(z) \{{z-Zk) nmk + Y. Y[(z-Zk)- {n+1)mk +r N (z), 

n-0 k=l n=0 fe=l 

where q n ,m{z) is the polynomial of degree m — 1 

p TP (z — 7, ) mk rnj ~ 1 

(3.2) q n>m (z) = ^ =, -— ^ Onj,l[z - Zj) 

j=1 llk=l,k^\ Z 3 ~ Z k) [=0 

and the coefficients a n ,j.i of this polynomial are given by the Cauchy integral 

J_ f /H dw 

( ' nd ' 1 ~ 2m J Tl (w - WLii 1 * ~ z *) nmk ' 

Also, t n , m (z) is the polynomial of degree m — 1 

(3-4) t n , m (z) . £ llt hk % )mk E ~ *)', 

where the coefficients b n> j t i of this polynomial are given by the Cauchy integral 

(3.5) b ndll = -L / f[( W ~z k r^ \ fiw)d ™ . 

The remainder term Vn{z) is given by the Cauchy integrals 

/H dw TU*-*> ^ Nmk - 

k=l 



. . 1 f f(w)dw , 

r N{z)=— t \Ttp — 1 \lw ll'* - ^ 

2m J r (w - z)[[ F k=1 (w - Zk) Nmk f-J: 

(3.6) 



If 

— Y[( w - z k) 



Nm k f( W ) dw 1 



In these integrals, the contours of integration Ti and T2 are simple closed loops 
contained in ft \ Clo which encircle the points z\, z-z,...,z v in the counterclockwise 
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direction. Moreover, z is not inside T 2 , whereas Ti encircles T 2 and the point z 
(see Figure 4 (a))- 

The expansion is convergent for z inside the "lemniscate annulus" (see Figure 
5) 

p 

(3.7) A v = {z e n \ n 0l r 2 < IJ \z - z k \ m " < n}, 

k=l 

where 

(3.8) ri = lrd we€ \ Q |jj \w - z k \ m > j , r 2 = Sup tt£fio j J[ \w - z k \ m " j . 
Proof. By Cauchy's theorem, 

(3.9) f(*) = — I f ^ dw 1 f /(w) dw , 

where Ti and T 2 are the contours defined above. First we assume that the m points 
of the set S are all distinct and later we will let the first mi points to coalesce to 
z\, the second m 2 points to coalesce to z 2 and so on. We substitute (|2.12|) - (|2.13|) 
into the first integral above and 

/o 1f y, 1 H m {w,z;zi_,...,z m ) 1 _ Y\£ =1 {w-Zk) 

w-z llfe=i(2-^fc) 1-u llfe=l(^-^j 

where H m (w, z;Z\,..., z m ) is defined in l|2.2|) . into the second one. Now we intro- 
duce the expansion i|2.1 4fl of the factor (1 — it) -1 in both integrals in i|3.9[l . Using 
(12. 3() and after straightforward calculations we obtain 

N — l m N—l m 

(3.11) }{z) = J2 1n,m(z) H(Z- Z k ) n + t n<m (z) JJ(« + 
n=0 fc=l n=0 fe=l 

where q nyT n(z) is given by formulas l|2.7|l and (|2.8|l replacing the contour C by IV 

Also, 

(3.12) 



=g^ nr=w . (z ._ Zfc) . ^ = ^y ra n(«) 



_ 1 f(w)dw 



and 



(3.13) 



2m J r (w - z) \\ k=x {w - z k ) N AX 



1 /■ m 



N f(w)dw 1 

N ' 



2ni Jr 2k l J{ " w ' z Uk=i( z ~ z k) 
Now we write 



(3.14) t n>m (z) = — Y[(w-z k ) n+1 f(w)dw^2 

JV9. 1. 1 1 



Uk=i,k&( z - z k) 



and repeat the steps following 1)2. 32|) in the proof of Theorem 12 . 51 for q n , m {z) and 

tn , m ( ^ ) • 
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FIGURE 4. The casep = 3. (a) Contours Ti and T2 in the integrals 
(|3.3[1 . I|3.5fl and (|3.6|> . (b) For z £ 4 p , we can take a contour T2 in 
51 located between fio and A p and a contour Fi in fl such that A p is 
inside this contour. Therefore, Yl k=1 l z — z fe| m,c < IIaUi — -z/c| mfc 

v ™ g ri and nLi l w - ^| m * < n*=i I* - ^l" lfc Vfer 2 . 




(r!,^) (r!,rf) (r?,r£) 

Figure 5. Shape of the "lemniscate annulus" A p for p = 3. It 
depends on the relative size of the parameters r\ and r2 defined in 
(|3.8|l . The different forms are labeled by (ri,T2) with r\ > r%. In 
these pictures \z 2 — z 3 \ < \z\ — z 3 \, \zi — z 2 \ and r% > r\> r\. 



For any z G A p we can take simple closed loops Ti and T2 in O \ Oo such that 
(see Figure 4 (b)) 

p p 

(345) Y[\z-z k \ m « < Y[\w-z k \ m " VweT 1 

k=l k=l 

and 

p p 

(346) Y[\z-z k \ m " >Y[\w-z k \ m " v w er 2 . 

k=l k=l 

On these contours |/(w)| is bounded by some constant C: |/(u>)| < C. Introducing 
these bounds in (|3.t)|) we see that liniAr^oo tm{z) = and the proof follows. □ 
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If the only singularities of f(z) inside f2o are just poles at z\, Z2,---,z p , then 
alternative formulas of (|3.3|) and (|3.5|l for computing the coefficients of the above 
multi-point Laurent expansion is given in the following proposition. 

Proposition 3.2. Suppose that gk(z) = (z — Zj,) Pk f(z), k = 1,2, ... ,p are analytic 
functions in Q for certain p\, pi,---,Pk G N. Then, for n = 0, 1,2, . . ., coefficients 
a n,j,i and &n,i,i vn expansion are also given by the formulas: 



(3.17) 



and 
(3.18) 



nmk+Pk-l 



gk{w) 



jynmj+pj+l 



(w-z^UU^ 



IE= Wi («-«.) T - 



^ ' £)Pk-(n+l)m k -l 



9k{w) 



(W - Zj) l + 



t n («-*.) (n+i)r - 



+ 



9j{w) n (w-2 S 



(n+l)r 



Proof. We deform both contours Ti and T2 of equations (|3.3(l and 13.5(1 . respectively, 
to any contour of the form CiUC^U- • -UC P contained in Q, where Cfc, fc = l,2,...,p 
is a simple closed loop which encircles the point z^ in the counterclockwise direction 
and does not contain the point Zj j = 1,2, ... ,p, j ^ k inside (see Figure 3 (b)). 
Then, 

(3.19) 



1 p 

— y 



9k{w) 



dw 



gj(w) dw 



and 



(3.20) 



IlLi,^ fc (^-^) ( " +1)m - 9k(w)dw 



J n,j,l 



2m 



E 



hi n 



(to — z fe Vfc-("+ 1 ) m fc 
gj(w)dw 



f w _ z _.)Pj-(n+l)mj+i+l 



□ 



From here, equations (|3.17(l and (|3.18|l follow 

Remark 3.3. Let z be a real or complex variable. Suppose that gk{z) = (z — 
Zk) pk f(z) is pfc — 1— times diffcrcntiable at every Zk in S for some pk £ iV, /c = 
1, 2, . . . ,p. Define 

Af p 



(3.21) 



g(z)=f( Z )-J2tn,m(z)l[( 



z - z k 



-(n+l)m fc 
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where M = |_Max{(pi - l)/ mi, { pi - l )/m 2 , ...,{p p - l)/m p }\ and t n ^ m {z) is the 
polynomial defined in l|3.4(l and (|3.18(l . Then, the thesis of Proposition 12.91 holds 
for f{z) replaced by g(z). Moreover, if Ofc=i( z — z k) Pk f{z) is an analytic function 
in f2, then the thesis of Theorem 12.51 applies to g(z). 

4. Multi-point Taylor-Laurent expansions 

For multi-point expansions we have the possibility (that we do not have in the 
standard theory) of expanding in Taylor series in some points and in Laurent series 
in other points. 

Theorem 4.1. Let SIq and Q be closed and open sets, respectively, of the complex 
plane, and £!o C C (D. Let f{z) be an analytic function on f2\f2o> Z\, Z2, ■ ■ ■ , z q £ 
fl\flo and z q+ i, z q+ 2, ■ ■ ■ , z p G CIq (l points are in andp — q points are in £1$). 

Write s = mi + mi + • ■ • + m q . Then, for z £ Cl\ flo, f{z) has the Taylor- Laurent 
expansion 



n=0 k=l n=0 nk=q+l\ Z Z k) 



(4.1) 

Vf q (7 - 7, \(n+l)m k 

■ ni, + i(^%) (n+lK 

where q n ,m{z) is the polynomial of degree m — 1 

p n p - ( z - z k) mk mj ~ 1 

(4.2) q n>m {z) = T[V~ M (z--zu )^ E a ^ z - z i? 

? '=1 iik=l,k^j\ 3 k) ;=0 



and the coefficients a n j.i of this polynomial are given by the Cauchy integral 
_ 1 f / M d w 

andtn} m {z) 

respectively, 



Also, tn}m(z) andtn)m(z) are the following polynomials of degrees s — 1 andm—s — 1 



(4.4) ew ^-H wiz *U e h ^ ~ 

where the coefficients b n jj °f this polynomial are given by the Cauchy integral 
(4.5, KtJ*— ' D^t!'-^ «"»*» 



2iriJ r! IlLiC 10 - z k) nm " (w-Zj)' +1 ' 

(4-6) ew - E E - *)'. 

where the coefficients c n ji of this polynomial are given by the Cauchy integral 

U7) _ 1 f UL q+ A^~ z k) in+1)mk fHdw 

(4. I ) C n jJ 



2ir!i r! Ul=i(™ ~ z k) (n+1)mk {w-Zj) l + 1 ' 
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The remainder term tn(z) is given by the Cauchy integrals 



tn{z) 



(4.8) 



f(w) dw 



, jVrrifc 



2m J Fl (w - z) lH =1 (w - z k ) N ™» **» 



Wr 2 nLi(^-^) Wmfc itc=«+i (*-**) 



Ni. 



In these integrals, the contours of integration T\ and T 2 are simple closed loops 
contained in Q \ Qq which encircle f^o in the counterclockwise direction. Moreover, 
the points z and Z\, z 2 ,...,z q are not inside T 2 , whereas I\ encircles T 2 and the 
points z and Z\, z 2 ,...,z q (see Figure 6 (a)). 
The expansion is convergent in the region (Figure 7) 



D q , p = \ z e Q \ fi , [J |(2 - z k )\ m * < n, 



(4.9) 



fe=i 

fc=i 



z k )\ mk <r 2 J] K*"**) 

k=q+l 



where n = Inf we e\n {Hl=i \( w ~ z k)\ mk } and 



r 2 = Inf weno {]JUi \( w ~ z k)\ mk U' k=q +i H"' ~ -/. 
Proof. By Cauchy's theorem, 



(4.10) 



/(*) - 7T-- 

2m 



f{w) dw 



w — z 



1 

2ni 



f{w) dw 



w 



where T± and T 2 are the contours defined above. 

First we assume that the m points of the set S are all distinct. Later we will 
let the first m\ points coalesce to z\, the second m 2 points to z 2 , and so on. We 
substitute (|2.12() - (|2.13() into the first integral above and 



(4.11) 
where 
(4.12) 
and 

(4.13) 



1 



F m (w j z, Z\ , . . . , 077J.) 



1 



]lfc=l( w - z k) IL= S +1 & -Z k )l-U 

Uk=i( z - z k)UkL s +i( w - Z k) 

' u s k=i( w - z k) nr= s +i( z - z k) 



Fmi.1V , Z, Z\, . . . , Z m ) 



W 



Y[(w-z k ) Y[ (z-Zk)- 

_fc=l k=s+l 
s m 
y^iz-Zk) ]J[ i w ~ z k) 



k=l 



k=s+l 



into the second one. Next we introduce the expansion Ij2.14|l of the factor (1 — u) 1 
in both integrals in l|4.1()JI . We observe that F m (w, z; z\, . . . , z m ) may be written as 
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F m (w,z;zi,...,z m ) =H s (w,z;zi,...,z p ) JJ_ (z - z k )~ 
(4.14) fe=s+1 

[w, z; z p+ i, . . . , z 

III 

)\{{z~z k ), 

where H m (w, z; Z\, . . . , z m ) is defined in l|2.2ll . Using this decomposition, equation 
(I2.3|> . and after straightforward calculations we obtain 

N-l m N-l t-js / _ \ n 

f( Z ) = e qn , m (z) \{{z- z k r + E CW irH 1 , n + 

1C , n=0 fe=l n=0 llfc=s+l^ Z ^ 

(4-15) JV _ 1 

E,(2) IL^i^ ~ z k) n , 
b n,mK*) j-jm + l ( z _ Zfe )"+1 W ^ ^' 

where q n:in (z) is given by formulas i|2.7|) and (|2.8(1 replacing the contour C by Ti. 
Also, 

(4.16) =-E^ £ Ufc ^ ( r^ 

J = 1 llfe=l,fc^il Z J Z k) 

with 

(4.17) 1 f nr^ +1 (— ^)-/h^ 



27r *Vr 2 IlLiH-^)" 

m 

(4-18) 4 2 LM = E 



rifc^a+i.fc^^ z fc) 

n J nm 7 \" 

j=s+l llfc=s+l,fc^jl Z j — z *v 



C 



with 

(4.19) 

and 

(4.20) 



i r uT= s+ i^-^r +i f(w)d W 



71,1 27T»7 ra nu(^-^)" +1 



^)=^/ 7 — — wfl^-^r- 

2m J Tl (w - z) \A k=l {w ~z k ) N 

w *V„,^ ,7,., TT S ^ » \N 



Now we write 
(4.21) 



2 "ir 2 nLi( w - z fe) n ^ ( w - z j)Ul=i,kjtj( z 3 - z k) 

, (2) (z) = ±f nL +1 (^-^)" +1 v dw 



J = S + 1 

(1) 



and repeat the steps following l|2.32|l in Theorem 12.51 for q nm (z), tn' m (z) and 
t {2) (z) 
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(a) 



(b) 



FIGURE 6. The case q = 2, p = 3. (a) Contours Ti and T 2 in the 
integrals (|4~!)|) . (f4~7)) and (02J. (b) For z £ we can 

take a contour r 2 located between D,q and -Dq.p and a contour Ti 
in f2 with D qp inside this contour. Therefore, [J?j \z — Zk\ mk < 

IlLi \w-z k \ m * vw e r a andnLi k-^r* ]lL ? +i l«>-«*r* < 



m=ii^-^r n 



k =q+l 



z k \ m " v w e r 2 . 



For any 2: verifying 1)4. 9|l . we can take simple closed loops Ti and r 2 in S7 \ fio 
such that (see Figure 6 (b)) 

f[\z-z k \ m * < f[\w-z k \ m * VweTi 



(4.22) 
and 
(4.23) 



n 



fc=i 



w - z k \ 



fc=i 



n 

k=q+l 



Z k \ 



g 

>n 

k=l 



z k 



Yl \ w - z k\ 



k=q+l 



V w 6 r 2 . On these contours |/(w)| is bounded by some constant C: |/(w)| < C. 
Introducing these bounds in (|4.8|) we see that lim^^oo r^(z) = and the proof 
follows. □ 

If the only singularities of f(z) inside flo are just poles at z q+ i, z 9+2 ,...,z p , then 
alternative formulas of l|4.3l) . (|4.5|) and (|4.7|l for computing the coefficients of the 
above two-point Taylor-Laurent expansion is given in the following proposition. 

Proposition 4.2. Suppose that g k (z) = [z — z k ) Pk f(z) is an analytic function 
in for certain p k G N and k = q + 1, q + 2, . . . ,p. Define g k (w) = f(w) for 
k = 1,2,3, . . .q. Then the coefficients a n ,j,i, b n .j,i and c n jj in the expansion \4-l\ ) 
are also given by the formulas: 

(4.24) 



a ■ 1 — N jynm k -l 



fH 



k—q-\-l,k^j 



9k(w) 



■Qnmj+pj+l 



(^^■) I+l nii, s# (»-^)' 

9j{w) 
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Figure 7. The region D q p denned in Theorem 14.11 is given by 
D q .p — O p p| B q>p , where O p is the "lemniscate domain" of foci Zi, 
. . ., z p and parameter r%. Also, B q , p = {z G <D, Y\ 9 k= i \(z — Zk)\ mk < 
r 2 nfc= 9 +i l( z — z k)\ mk }- This pictures show the topologically dif- 
ferent forms of D q ^p depending on the relative value of r\ and r 2 
when q = 2 and p = 3. The pictures are labeled with (ri,r 2 )- In 
these pictures z\ \z\—z 2 \ < \z\ — z?\, \z 2 — z 3 \ and r 2 > r 2 > r 2 > r 2 . 



(4.25) 



°n,j,l 



k=q+l 



jjp k -nm k - 



' 9kWn p a=q+ltaftk (w-z s r m < 
(w-^) ,+i n2=i(«'-«.) nm * 



(4.26) 



p 



r . ,_ 7->Pfc-(n+l)m fc -l 
k—q+l^k^j 



£ ) Pj-(n+l)m j +; 



g fc Hnu+i,^(^-^) ( " +1)r 
(w-^o'^riLi^-^)^ 13 " 

9iWnUi^(»-^) ( " +1) ™ 



mu(™-o (n+i)r 



Proof. We deform both contours Ti in equation (|4.3() and the contour r 2 in equa- 
tions l|4.24|) and 1)4. 7[) into any contour of the form Ci U C2 U • • • U C p contained 
in Q, where C'k, k — 1,2, ...,p is a simple closed loop which encircles the point Zk 
in the counterclockwise direction with Zj not inside Ck, j = 1, 2, . . . ,p, j ^ fc (see 
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Figure 3 (c)). Then, 
(4.27) 



1 f 9k{w) dw 

2m JzU, Jc k («> - z 3 ) l+1 T\ P s =i,sM w - Zs) nm » (w - z k ) nm * + r* ' 
gj (w) dw 



<i ]lLl, s #i( w - Z s ) nm ° (w - z^nmj+pj+l+l ' 
(A 08"\ h 1 f ITLg+M^fcO " Z s) nms 9k{w)dw 

(4.28) b n<jtl - — ^ 



2wi J Ck (w - ZJY+ 1 Ut=i(w ~ z s ) nm ° (w - z k )Px- nm * 



(4.29) 

071,3,1 kM.^ 2m ' c » {w z > y+1 n - i(w - ^) (n+i)ms ( w - zky*-^™* 

From here, equations 14.2411 . (|4.25|l and (I4.26[) follow. □ 

Remark 4.3. Let z be a real or complex variable and suppose that (z — Zk) pk f(z) 

is pk — 1— times differentiable at Zk for certain € N. Define 

(4.30) 

j£ n 9 (7-7, \ nm >> FT 9 (7 - 7, \(n+l)m h 

g{z)^f{z)~Yt£Uz) Jil k = l[ , "\ Yt%U*) ™, u+n > 

n=0 \.\.k=q+l\ Z Z k ) n=Q llk=q+l[Z Z k )y 

where M = LMax{(p g+1 -l)/m g+ i, (p q+ 2 ~ l)/m q+2 , ... , (p p -l)/m p }\ and t$ m (z) 
and tn} m {z) are the polynomials defined in 1)4.4(1 . I4.6jl . I|4.25(l and l|4.26j) . Then, 
the thesis of Proposition 12.91 holds for f(z) replaced by g(z). 

Moreover, if nfe= g +i( z ~ z k) Pk f(z) is an analytic function in f2, then the thesis 
of Theorem 12 . 51 applies to g(z). 

5. Discussion and concluding remarks 

In an earlier paper |H] we have discussed the theory of two-point Taylor expan- 
sions, two-point Laurent expansions and two-point Taylor-Laurent expansions. In 
the present paper we have generalized these two-point cases to multi-point cases. 
We have given details on the regions of convergence and on representations of the 
coefficients and the remainders of the expansions in terms of Cauchy-type integrals. 

Multi-point Taylor expansions are related with topics from interpolation theory, 
in particular with the Newton interpolation theory with applications in numerical 
analysis. For example, applications can be found in initial and boundary value 
problems in connection with ordinary differential equations and in numerical quad- 
rature of integrals. 

From the point of view of interpolation theory detailed information on multi- 
point expansions can be found in Chapters 3 and 8. The theory of several-point 
Taylor expansions is discussed in Chapter 3 of [H], although in a setting that is 
different from our approach. Our approach gives explicit Cauchy-type integrals 
of coefficients and remainders which cannot be found in Walsh's approach. In 



MULTI-POINT TAYLOR EXPANSIONS OF ANALYTIC FUNCTIONS 



1!) 



particular, we cannot find explicit formulas for the polynomials q n ,m(z) of formula 
(15) as we have in (16)-(17). Knowledge of these explicit formulas is necessary to 
construct asymptotic expansions of integrals with several saddle points. 

In addition to this, our Laurent and Taylor-Laurent expansions are new. They 
have a formal similarity with the rational approximations of Chapter 8 of [8]: they 
involve negative powers of z. But they are completely different. The rational 
approximations, in particular the Pade-type approximations P n (z)/Q m (z) are of 
interpolatory type. These are generalizations of the Taylor polynomial at several 
points: a quotient of polynomials instead of a polynomial. However, our expansions 
(21) or (32) have a different form and a different approximation property: they 
approach not only at regular points like Pade-type approximations but also at 
singular points of f(z). And of course, the regions and convergence properties in 
[Hj are different from ours. 

Apart from applying the present results in problems from interpolation theory, 
in particular in problems from numerical analysis, we expect to find applications in 
asymptotic analysis of integrals, which application area is our main motivation; see 
[5] . In that paper certain orthogonal polynomials have been considered and we have 
given new convergent expansions that also have an asymptotic property for large 
values of a parameter (the degree n of the polynomials). Orthogonal polynomials 
and special functions can be studied when the variable and several parameters are 
large. In that case more than one or two so-called critical points occur that may give 
the main contributions to the integral, and expansions of analytic functions at these 
points gives again the possibility of constructing new convergent expansions with 
an asymptotic property. This method avoids the complicated conformal mapping 
of the phase function of the integral into a standard form (say a cubic or higher 
polynomial). In addition, when the critical points are multiple poles, Laurent- 
type expansions may be considered. A few application areas are mentioned in the 
Introduction, see the integral in (|1.4Jl , which we expect to approximate in terms of 
Airy functions and the Pearcey integral l|1.5J) and its derivative with respect to x 
and y. 
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